Shot Noise Probing of Magnetic Ordering in Zigzag Graphene Nanoribbons 
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The nonequilibrium time-dependent fluctuations of charge current have recently emerged as a 
sensitive experimental tool to probe ballistic transport through evanescent wave functions introduced 
into clean wide and short graphene strips by the attached metallic electrodes. We demonstrate that 
such "pseudo-diffusive" shot noise can be substantially modifled in zigzag graphene nanoribbon 
(ZGNR) due to the topology of its edges responsible for localized states that facilitate ferromagnetic 
ordering along the edge when Coulomb interaction is taken into account. Thus, the shot noise 
enhancement of unpolarized, and even more sensitively of spin-polarized, charge currents injected 
into ZGNR will act as an all-electncal and edge-sensitive probe of such low-dimensional magnetism. 

PACS numbers: 73.50.Td, 75.75.-|-a, 73.63.-b, 81.05. Uw 



The advent of graphene [ij — first truly two- 
dimensional crystal whose carbon atoms form a 
honeycomb lattice — has reinvigorated exploration of 
low-dimensional quantum transport phenomena. The 
high-mobility graphene flakes are far more accessible 
to different experimental and fabrication techniques 
than a traditional two-dimensional electron gas (2DEG). 
Furthermore, the chiral massless Dirac fermions carrying 
the current in graphene are quite different from the 
usual quasiparticles in 2DEGs. Thus, many standard 
mesoscopic transport experiments have been reexam- 
ined in this new 2D setting to unveil their anomalous 
versions [H due to relativistic-like character of low 
energy electronic excitations. In particular, very re- 
cent experiments on the shot noise in ballistic 2] and 
disordered Q graphene strips elucidate fundamental 
conduction properties of Dirac fermions that cannot be 
extracted from conventional time-averaged current /. 

The shot noise denotes time-dependent current fluctu- 
ations, driven by the nonequilibrium state, which persist 
down to zero temperature and originate from the gran- 
ularity of charge The zero- frequency noise power S 
and related Fano factor F = S/2eI can probe the effects 
of disorder, carrier statistics, and interactions in sam- 
ples smaller than the electron-phonon inelastic scattering 
length 0, |3| . The Poissonian limit _F = 1 characterizes 
transport governed by uncorrelated stochastic processes 
(as encountered in tunnel junctions). More intriguing 
sub-Poissonian Fano factors are found in, e.g., noninter- 
acting diffusive conductors where F = 1/3 is determined 
by the interplay of quantum stochasticity due to impurity 
backscattering and the Fermi statistics 

Surprisingly enough, it has been predicted Q that the 
same F — 1/3 should be measured in clean wide and 
short graphene strips {W/L > 4 where W is the width 
and L is the length of the strip). This is in sharp con- 
trast to ballistic transport (with perfect transmission) 
in 2DEG channels where shot noise is absent (F = 0) 
due to completely correlated propagation of electrons 
by the Pauli principle. Such "pseudo-diffusive" Q con- 
duction at the Dirac point of graphene can be traced 



to metal- induced gap states [6|], familiar from metal- 
semiconductor junctions where they provide direct tun- 
neling between evanescent scattering states of external 
leads populating the gap of a short sample. However, 
in graphene, viewed as a gapless semiconductor where 
valence and conduction bands touch at the Dirac point 
{Ep = 0) [H , the corresponding evanescent states [3] can 
penetrate much longer distance carrying current whose 
F = 1/3 is accidentally 0| the same as in the diffusive 
metallic wires. In fact, F = 1/3 has been observed in 
recent experiments on large aspect ratio {W/L ~ 24, 
L ^ 200 nm) ballistic two-terminal graphene devices 

We recall that F — 1/3 for two-terminal disordered 
conductors is universal — it does not depend on the im- 
purity arrangement, band structure, and the shape of 
the sample [4] . Much less is know about the universality 
of the "pseudo-diffusive" shot noise in ballistic graphene 
strips where typically GNRs with armchair edges [a, |3, 0| 
or infinite gated ZGNR setups have been the subject 
of theoretical analysis. The topology of the zigzag edge 
is rather special generating peculiar edge-localized quan- 
tum states ll|. Moreover, their partially flat (within 
one-third of ID Brillouin zone) subband generates a large 
peak in the density of states at the Fermi energy Ep = 0. 
This instability [l3| is most likely resolved [13| through 
magnetic ordering around the zigzag edge when electron- 
electron interactions (even infinitesimally small [lH) are 
"turned on." Such carbon-based magnetism involving 
s-p orbitals was conjectured for an infinite ZGNR us- 
ing the Hubbard model [flj , and confirmed through nu- 
merous recent density functional theory (DFT) calcula- 
tions 13, 14,[lB| and related proposals for spintronic de- 
vices |16l. Il7|. 

Here we address these unresolved issues by predicting 
two experimentally testable effects in clean two-terminal 
mctal-ZGNR-metal devices that strongly intertwine their 
magnetic correlations with the electrical shot noise: («) 
finite length ZGNR attached to metallic electrodes will 
also develop magnetic ordering around the edges which, 
however, decays in the vicinity of metallic contacts and 
requires finite strength of Coulomb interaction (Fig. [2]); 
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(n) the magnetic moment per carbon atom in the mid- 
dle of ferromagneticaUy ordered zigzag edge is in one-to- 
one correspondence with the enhanced shot noise Fano 
factor F > 1/3, so that all- electrical measurements can 
be used as edge-sensitive technique to probe unusual s-p 
magnetism in wide ZGNRs at low temperatures (Fig Bt. 

We employ the single 7r-orbital Hubbard model to 
obtain magnetic ordering within ZGNR: 
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which is defined on a finite-size honeycomb lattice with 
zigzag edges and the lattice constant a — 2.46 A. Here Cj 
(ci) creates (annihilates) an electron in the 7r-orbital lo- 
cated at site i = {ix,iy) and 7 is the nearest neighbor hop- 
ping. The width of the A^^-ZGNR lattice is measured us- 
ing the number of zigzag longitudinal chains [HI, [3] ■ 
We use the number of atoms comprising a zigzag 
chain to measure its length [in the units of a, the average 
width of ZGNR is = a^/^{N^ - l)/2 and its length is 
L = {N^- l)a/2]. The metallic leads attached to ZGNR 
are modeled by the tight-binding Hamiltonian {U = 0) 
on the semi-infinite square lattice [7 J8[ , which allows us 
to capture different contact effects 7| introduced in ex- 
perimental circuits [l|, 0, 0| by ultimate electrodes being 
metals rather than all-graphitic structures. 

By setting U — both in the electrodes and in the 
ZGNR central region, the Fano factor in Fig. [TJa) is ob- 
tained using the celebrated scattering theory formula [3| , 
F = EnTnil - r„)/^„r„. Here T„ are the eigenval- 
ues of ttt and t is the transmission matrix of a phase- 
coherent device. In the basis of eigenchannels that di- 
agonalize tt^, a mesoscopic device can be viewed as a 
parallel circuit of independent one-dimensional conduc- 
tors. Figure [TJb) shows that Fano factor is substantially 
affected by the type of the square lattice leads 0, even 
in samples with large W/L. They effectively introduce 
disorder at the lead/ZGNR interface, thereby mixing the 
transverse propagating modes. That is, t acquires non- 
zero off-diagonal elements even in clean devices, as shown 
in Fig. [UJb). The absence of mode mixing in trans- 
port through armchair GNR, attached to highly doped 
graphene 5] or square lattice metallic leads JTI, Jsj, l9| , is 
intimately connected to its F = 1/3 value [1, HIJ- Nev- 
ertheless, for metallic leads in Fig. [Ija), whose lattice 
matches f?\ the honeycomb lattice (square lattice spacing 
is the same as the carbon-carbon distance) while allowing 
propagating modes in the leads to efficiently couple via 
evanescent ones in ZGNR in Fig. H^b), we find Ohmic- 
like conductance G = 2|! 7;^ _ w/L in Fig.^a) and 
the Fano factor ~ 0.3 - 0.4 for W/L = 5 in Fig. [T^b). 

We select 212-ZGNR of length = 75, whose di- 
mensions ensure F ~ 1/3 in Fig. [Ijb), for the analysis 
of interaction U > driven magnetism and its effects 
on the shot noise. In infinite ZGNRs, edge states mani- 
fest as two peaks in the local density of states along the 
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FIG. 1: (Color online) (a) The Fano factor and conductance 
of nonmagnetic {U — 0) 212-ZGNR attached to two metal- 
lic electrodes modeled as square lattices, (b) Shot noise in 
ZGNRs of variable width W and length L, but with their 
ratio fixed at W/L = 5, for two diflferent types of metallic 
electrodes whose square lattice is "matched" (solid line) or 
"unmatched" (dashed line) to the honeycomb lattice. The 
electrons are injected at the Fermi energy Ef = 10~^7. 



edge [18|. The overlap of states from two edges yields 
bonding and anti-bonding states enabling a single con- 
ducting channel close to the Dirac point Ep — with 
highly unusual transport properties fl8]. The evanes- 
cent state enabled Ohmic-like transport G oc W/L in 
finite ZGNR devices is expected at energies Ep < AE 
of the single open conducting channel [10]. Within non- 
interacting {U = 0) 212-ZGNR this channel is the only 
open one when Ep of the injected electrons is Ep < 
AE — 0.015i while 317 channels are used for in- 
jection from the square lattice lead in Fig. [TJa) at half 
filling. The Fano factor and conductance of this ZGNR 
are determined by few non-negligible transmission eigen- 
values, out of which one is close to unity corresponding 
to transport through the edge states induced channel. 

To capture the magnetic ordering in wide ZGNR two- 
terminal device, we use the standard mean-field decou- 



pling scheme [l2] which yields an effective ll|,ll7| single- 
particle Hamiltonian 
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The inhomogeneous exchange potential Ai — ~U{(hi^) — 
(nil)) determines the spin splitting, where a^^ are ele- 
ments of the Pauli matrix. The local magnetic ordering 
is quantified by the z-component of spin (S?) = {{fii^) — 
{nil))/2 and the local magnetization Wj = gfisiSi). The 
average electron density in equilibrium 
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Tii^) = dElm{ia\G''{E)\ia)fiE - Ep)l3a) 



G^{E) — [E — Hmfa — '^1 — '^2] ^, 



(3b) 



is computed from the diagonal matrix elements of the 
retarded Green operator G^{E) for the open system 
ZGNR-I- leads [l^l- The retarded self-energies operators 
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El, S2 introduced by the interaction with the leads de- 
termine escape rates of electrons into the electrodes. 

Equations ([2]) and ([3]) are solved self-consistently with 
charge convergence accelerated by using small non-zero 
temperature in the Fermi function f{E — Ep) describ- 
ing the electrodes [31 . The result for the local magne- 
tization plotted in Fig. [Ifa) shows usual staggered pat- 
tern with salient feature being ferromagnetically ordered 
spins along each edge and opposite spin directions be- 
tween the edges as the ground state [HI, O, [3] . How- 
ever, Fig. m^a) highlights two major differences between 
finite ZGNRs sandwiched between two metallic leads and 
infinite ZGNR: («) the magnetic moment per edge car- 
bon atom decays upon approaching the metallic lead, so 
that no magnetic ordering appears in ZGNR shorter than 
~ 6; [ii) while infinitesimally small U causes ferri- 
magnetic ordering on each sublattice around the edges 
of an infinite ZGNR [llj, in two- terminal ZGNR devices 



[7/7 > 0.2 is required to get sizable modgc on outermost 
carbon atoms, as shown in the inset of Fig. [ll^a). 

Although DFT calculations go beyond only on-site 
Coulomb interaction and nearest-neighbor hopping of 
-ffiviFA, they can be mapped [isl . [l9j to electronic and 
magnetic structure at half-filling (we use Ep — 10^^7) 
obtained from simpler Eq. ([2]) . The replacement of i?DFT 
with minimal-basis-set Hmfa allows us to treat systems 
composed of ~ 10^ carbon atoms. The values for U and 
7 extracted via this mapping depend slightly on the type 
of approximation scheme for exchange-correlation den- 
sity functional [H. We find that U/j = 1 in Fig. d^a) 
reproduces the ab initio result [Tsj l for the magnetic mo- 
ment ?7i"J[g^ « 0.28 fiB per carbon atom in the middle 
of a zigzag edge, as well as the energy gap A" [19] of 
an infinitely long ZGNR opening around Ep — due 
to the staggered sublattice potential [ij. The gap A° 
is experimentally visible [20| in very narrow ribbons and 
vanishes [T^ within the room-temperature thermal en- 
ergy window when the width of ZGNR reaches ~ 80 nm. 

The self-consistently computed effective Hamiltonian 
Eq. ^ is unaffected by small voltage bias [l^ and can 
be used as the basis for linear response transport calcula- 
tions within single-particle formalisms 17[ . While recent 
studies have employed such Hamiltonian to obtain the 
spin-resolved conductances of two-terminal ZGNR de- 
vices [l3] , here we utilize the scattering approach to quan- 
tum transport [3| to compute the spin-dependent shot 
noise as quantified by correlators between spin-resolved 
charge currents /J and I2 in the drain lead 2: 
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si^'it')si^it)). (4) 



Here /^(t) is the quantum- mechanical operator of spin- 
resolved charge current of spin-a electrons in lead 2. 
The current-fluctuation operator at time t in lead 2 is 
Si^ (t) = (t) ~ (i^ (t)) and (. . .) denotes both quantum- 
mechanical and statistical averaging f3| over the states 






[7/7 = 


= 1 ' 
















i 






FIG. 2: (Color online) (a) The spatial profile of the local 
magnetization mi (pointing out-of-plane) within 212-ZGNR 
of length Na ~ 75 {W/L = 5) induced by on-site Coulomb 
interaction U/'y = 1. The metallic leads are modeled as semi- 
infinite square lattice attached at the top and bottom arm- 
chair interface, (b) Amplitude of the transmission matrix 
elements \tqp\ for magnetically {U/'y = 1 panels) ordered 212- 
ZGNR from panel (a) and its nonmagnetic ([/ = panels) 
version studied in Fig. [T] using "lattice-unmatched" (left col- 
umn) or "lattice- matched" (right column) square lattice leads. 



in the macroscopic reservoirs to which ZGNR conduc- 
tor is attached via semi-infinite interaction-free metal- 
lic leads. The Fourier transform S22 (w) — J d{t — 
^1^^ ^-iu]{t-t ) gaa ^-j- _ ^/^ gives the spin-resolved noise 
power, and the total charge current noise is obtained 
from >5'22 ('-^) = So- '^'22^ (ijj) where we focus on the zero- 
frequency limit 5^2 — '5'22('^ ^ 0) of these expressions. 

The presence of intrinsic (due to exchange interaction, 
spin-orbit couplings, magnetic impurities, etc.) and ex- 
ternal magnetic fields in the sample is manifested in the 
shot noise dramatically when injected current is spin- 
polarized [U, [22, [2^ . This is due to the fact that any 
spin flip converts spin-t subsystem particle into a spin-| 
subsystem particle, where the two subsystems differ when 
spin degeneracy if lifted. Thus, the nonconservation of 
the number of particles in each subsystem generates ad- 
ditional source of current fluctuations. In particular, the 
shot noise of spin-polarized current injected from ferro- 
magnetic reservoirs into a ferromagnetic wire was pre- 
dicted to be a sensitive probe of its magnetic ordering 
and spin-dependent interactions (23j . 
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FIG. 3: (Color online) The Fano factor of magnetically or- 
dered two-terminal 212-ZGNR of length iV^ = 75 (VF/L = 5). 
Panel (a) is for unpolarized injection and two different types of 
metallic electrodes, whose square lattice is "matched" (solid 
line) or "unmatched" (dashed line) to the honeycomb lattice. 
Panel (b) plots total and spin-resolved shot noise for 100% 
polarized injection, where injected spins are aligned with the 
transport direction (the x-axis) and both spins are detected in 
electrode 2. The Fano factor for partially polarized injected 
current in the same setup as (b) is shown in (c). 



The computation of S22 for the general situation, when 
injected spins are characterized by a density matrix p^^, 
cannot be done using the standard basis of eigenchan- 
nels and the corresponding transmission eigenvalues r„ 
since in this approach the information about pf^ is lost. 
Instead, one has to use the lengthy formulas, derived in 
Ref. [2IJ in terms of the spin-resolved transmission matrix 



^21 
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— IniEj^' of the device, 
to keep track of pf^^ = (1 + Pin • o-)/2 and its spin po- 
larization vector Pin. Although no true long-range order 
is possible in one dimension (the spin correlation length 
along the zig zag edge decays to ~ 1 nm at room tem- 
perature [l5|), the observation of the shot noise requires 
low temperature 0] where magnetic correlations should 
be highly visible. Therefore, we use T = scattering for- 
mulas [21I [i^ l for S22 ■ The electrode 2 is paramagnetic 
so that it collects both spin species. 

The Fano factor F(Pin) = S^/2el2 for unpolarized 
I Pin I = current injection from electrode 1 into ZGNR is 
shown in Fig.[3ja). In the case of "lattice-matched" leads, 
the sensitivity of transport around Ep = to edge poten- 
tial [and reflection from it, Fig.[2jb)] drives slight increase 
[F ~ 0.37 at [7/7 = 1) of the Fano factor above the ref- 
erence value F ~ 0.35 at C/ = 0. The enhancement — 
F ~ 0.95 at U/^ = 1 — is much more dramatic for 
"lattice-unmatched" leads, where Fig. El^b) shows that 
even at C/ = propagating modes in two leads are de- 
coupled due to poor matching (the conductance does not 
scale as ~ W/ L) to most slowly Q decaying evanescent 
modes. At C/ > 0, a small gap (A^ ~ O.OO97 at U/"f = 1) 
opens in 212-ZGNR around Ep — which drastically 
reduces the transmission matrix elements [Fig. [2lb)], 
thereby, converting transport into conventional tunnel- 
ing with Fano factor close to F = 1. 



cross-correlators are zero S. 
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0. The shot noise of 



spin-polarized current becomes non-trivial when injected 
spins are non-coUinear with mi, as is the case of the x- 
axis (direction of transport) polarized spins in Fig. [3]^b). 
Since they are not the eigenstates of the local magnetic 
field along the z-axis, they are forced into precession . 
This enhances auto-correlation noise and generates 
non-zero cross-correlators 5*22 '"^ < 0. FigurelSfc) suggests 
that probing of ZGNR edge magnetic ordering via spin- 
dependent shot noise can be efficient even with partially 
polarized < [Pinj < 1 injected current. 

Conclusions — Our first principal result describes mag- 
netic ordering in wide and short ZGNRs where magne- 
tization persists at the edge (with no interedge mag- 
netic order characterizing narrow ZGNRs) even as its 
width goes to infinity. The edge magnetization dimin- 
ishes in the longitudinal direction upon approaching the 
ZGNR|metallic-electrode interface. This picture com- 
plements recent DFT calculations 13, 14, [l^ on nar- 
row infinite ideal ZGNRs. Since in the limit of wide 
ZGNR carbon-based s-p magnetism becomes a pure edge 
effect — not detectable by bulk-sensitive techniques — we 
offer a recipe on how spin-dependent shot noise can 
be exploited as edge- sensitive and all- electrical probe of 
this phenomenon. That is, Fig. [3] suggests an experi- 
ment that would: inject unpolarized current into ballistic 
ZGNR with ultrasmooth edges [2d| =J> measure possible 
enhancement of the shot noise above recently observed 
F = 1/3 inject in-plane polarized current while 
collecting both spins to observe further Fano factor en- 
hancement with increasing spin polarization [Pinj due to 
precession of transported spins in edge magnetic field. 
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